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ABSTRACT 

We examine plasma composition of jets in active galactic nuclei through the comparison of the 
total pressure (P) with partial pressures of electrons and protons in a cocoon. The total pressure 
is estimated from the analysis of an expanding cocoon dynamics. We determine the average kinetic 
energy per particle for several representative cases of particle energy distribution such as one- and 
two-temperature thermal plasmas and non-thermal electrons by evaluating the dissipation of total 
kinetic energy of the jet into the internal energy of cocoon plasma. The number density of the total 
electrons/positrons (n±) in the cocoon is constrained by using the particle supply from hot spots 
and the absence of thermal bremsstrahlung emission from radio lobes. By inserting P, n± and the 
particle energy of each population into the equation of state, the number density (n p ) and pressure 
(P p ) of protons in the cocoon can be constrained. Applying this method to Cygnus A, we find that 
(i) electron/positron (e ) pairs always dominate in terms of number density, but that (ii) either an 
"e ± -supported cocoon (i.e., P± > P p )" or "proton-supported one (i.e, P± < P p ) v is possible. 
Subject headings: galaxies: individual(Cygnus A) — magnetic fields — radiation mechanisms: non- 
thermal — radio continuum: galaxies — X-ray: galaxies 



1. INTRODUCTION 

Elucidating the formation mechanism of relativistic 
jets in active galactic nuclei (AGNs) is one of the greatest 
challenges of astrophysics in this century (e.g., Bland- 
ford and Znajek 1977; McKinney 2006; Komissarov et 
al. 2007; McKinney et al. 2012). Plasma composition 
of jets is a fundamental but difficult issue (Begelman et 
al. 1984 for review), because emission timescales of the 
bulk population such as low-energy electrons/positrons 
and protons are too long. To examine plasma composi- 
tion, discrete blobs in blazar jets have been utilized over 
the years. So far, three approaches have been pursued. 
The first is based on the synchrotron self-absorption limit 
combined with total kinetic powers of jets (Reynolds et 
al. 1996; Hirotani et al. 1999, 2000; Hirotani 2005). The 
literature indicates the existence of e ± pair plasma in M 
87, 3C 279 and 3C 345. The second is the constraint 
by the detection of circular polarization. Wardle et al. 
(1998) and Homan et al. (2009) examined the case of 3C 
279 and found that the minimum Lorentz factor of non- 
thermal electrons/positrons should be much larger than 
unity for electron-proton (hereafter e/p) content. They 
rather favored an alternative possibility of dominant e 
pair content with a small minimum Lorentz factor of non- 
thermal electrons/positrons (see, however, Ruszkowski 
and Begelman 2002). The third approach is the con- 
straint from the absence of bulk-Compton emission in 
flat spectrum radio quasars (Sikora and Madejski 2000; 
Ghiscllini & Tavccchio 2010) and it has been observation- 
ally tested for PKS 1510-089 and SWIFT J0746.3+2548 
(Kataoka et al. 2008; Watanabe et al. 2009). The same 
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approach has also been applied to the kiloparsec-scale 
knots in PKS 0637-752 (Georganopoulos et al. 2005; 
Uchiyama et al. 2005; Mehta et al. 2009). They claim 
that jets contain more e + e~ pairs than protons, but that 
jets are dynamically dominated by protons. However, 
it should be noted that the estimate of a total kinetic 
power Lj of each blob is difficult, because of the exis- 
tence of invisible components such as low-energy elec- 
trons/positrons and protons. Therefore, the assumption 
of constant Lj was made and the Lj arc inferred from non- 
thermal emissions. Since plasma composition is sensitive 
to Lj, a better estimate of Lj is essential. Regarding to 
the estimate of Lj, it is essential to take into account of 
the thermal component (e.g., Kino and Takahara 2008). 

Cocoons associated with Fanaroff-Rilcy I and II (FR I 
and FR II) radio galaxies are also known to be good tools 
for exploring plasma composition. In contrast to blobs 
in blazars, investigations using cocoon dynamics allow 
us to better estimate of energy injection into the cocoon. 
The total pressure P can be estimated with fewer uncer- 
tainties based on the dynamical interaction between jets 
and the intra-cluster medium (ICM) and P involves the 
contributions of invisible components (e.g., Rawlings and 
Sanders 1991; Fabian et al. 2002). For FR I radio galax- 
ies, many authors have discussed the ratio of P to that of 
non-thermal electrons (P^ T ) for various sources based on 
observed non-thermal emissions (e.g., Dunn et al. 2005; 
Croston et al. 2005; Rafferty et al. 2006; De Young 
2006; Birzan et al. 2008). First of all, we should empha- 
size that these studies indicate that the total pressure 
P tends to be larger than that of non-thermal electrons, 
i.e., P > P^ T (e.g., Birzan et al. 2008). This means 
that the finite pressure of low-energy electrons/positrons 
and/or protons is required in these sources. The derived 
P/pNT values in the 

previous work extend over a wide 
range from the order of unity to thousands (e.g., Birzan 
et al. 2008; Cavagnolo et al. 2010). For FR I sources, 



2 



Kino, Kawakatu, & Takahara 



however, an entrainmcnt process of surrounding medium 
via the jet boundary layer could work (e.g., De Young 
1993; Bicknell 1984; Rossi et al. 2008) and the process 
makes jets heavier. Therefore, jets in FR I sources could 
undergo severe proton loading during their propagations 
and this could cause the large scatter of P/P_. 

Instead, in this work, we focus on FR II radio galaxies 
(Fig. [1} from the viewpoint of the important advan- 
tage they represent. Contrary to FR I sources, we know 
from relativistic hydrodynamic simulations that no sig- 
nificant entrainment appears for FR II sources (Scheck 
et al. 2002; Mizuta et al. 2004). Therefore, a plasma 
composition test for FR II radio galaxies would allow 
us to give better constraints on plasma composition in 
AGN jets without an entrainment effect. Regarding an 
observational indication of a difference between total and 
non-thermal pressures in FR II radio galaxies, Ito et al. 
(2008) (hereafter 108) recently examined for the follow- 
ing sources (Cygnus A, 3C 223, 3C 284 and 3C 219). In 
108, they show that the energy density of total plasma is 
larger than the energy density of non-thermal electrons 
by the factor of 4-310 in the case of minimum-energy con- 
dition (e.g., Miley 1980; Kellermann and Pauliny-Toth 
1981). This implies that the minimum-energy condition 
is violated, particle energy is dominant, and low-energy 
electrons/positrons and/or protons (i.e., cosmic-rays) are 
required to explain the total P in these FR II sources. 

In §2, we describe the basic idea and assumptions of 
our method. In §3, we briefly explain the dynamical 
determination of the total pressure in the cocoon. In §4 
we express P as functions of the number density ratio 
of protons to electrons. In §5, we explain details of the 
plasma composition test. It is applied to Cygnus A in 
§6. Summary and discussions are given in §7. 

2. METHOD AND PROBLEM SETTING 

Here, we describe the basic idea and assumptions of our 
method. In this work, the number densities of protons 
(n p ), positrons (n + ) and electrons (ri-) arc related using 
the parameter 77 as follows: 

n p = rjn- 

n + = (l-r])n- (0 < 77 < 1), (1) 

where the latter relation is derived from the charge neu- 
trality condition. The case of rj = corresponds to 
pure e^ 1 plasma while r\ = 1 corresponds to the pure e/p 
plasma. We denote that n p = nJ + np T , n_ = +n^ T , 

n + = nT + /i+ T , and n± = ri- + n + where n± is the sum 
of the total number densities of electrons and positrons. 
Hereafter, superscripts T and NT represent thermal and 
non-thermal components, respectively. The distinction 
between thermal and non-thermal particles may not be 
trivial for relativistic plasmas. In this paper, we refer 
the thermal component to Maxwellian-like distribution 
characterized by the temperature, while we refer the non- 
thermal one to particles following a power-law distribu- 
tion characterized by the power-law index and minimum 
and maximum energies as detailed below. Since we fo- 
cus on relativistic plasmas in the present work, thermal 
component correspondingly has a relativistic tempera- 
ture. Hence one should be cautious since most of ob- 
servational papers refer the thermal component to non- 
rclativistic plasmas (e.g., Garrington and Conway 1991). 



The allocation of partial pressure of each plasma pop- 
ulation is the central concern of this paper. In general, 
P is decomposed to 

P = P- +P+ +P p + Pb 

= P T + Pj + F NT + P* T + Pj + F p NT + P B , (2) 

where Pi, P$, Pj, , P* T , P p NT , and P B are, 
the partial pressures of thermal (T) electrons, thermal 
positrons, thermal protons, non-thermal (NT) electrons, 
non-thermal positrons, non-thermal protons, and a mag- 
netic pressure respectively. We also define P± = P- + 
P + as the sum of the total pressures of electrons and 
positrons. Throughout this work, we do not include the 
magnetic pressure because it is sub-dominant in P. Isobe 
et al. (2005) summarize the energy density of energetic 
electrons as typically being 10 times larger than that of 
magnetic fields in various radio lobes (e.g., Isobe et al. 
2002; Tashiro et al. 1998, 2009; Hardcastle and Croston 
2010) and it also holds in Cygnus A (Yaji et al. 2010). 

2.1. Basic idea of the method 

The essence of our method is as follows. First, the total 
pressure in the cocoon (P) is determined through dynam- 
ical considerations following 108 where they obtained P 
via the comparison of the expanding cocoon model with 
radio observations. Second, average energy per one parti- 
cle in the cocoon is evaluated. It is essential that our for- 
mulation is based on the basic conservation laws of mass, 
momentum, and energy in the cocoon. Since it depends 
on coupling of protons to the electrons/positrons, we ex- 
amine several representative cases with different equa- 
tions of state. Third, n_ can be partially constrained 
by using the absence of thermal bremsstrahlung emis- 
sion from the cocoon and the supply rate of electrons 
from the hot spots. Finally, n p and P p can be obtained 
by inserting the obtained quantities into the equation of 
state (EOS). 

2.2. On particle distribution functions 

Since observational data at low frequencies below GHz 
are quite limited, it is hard to explore the properties of 
low-energy electrons (including positrons). Bearing this 
difficulty in mind, we pick up plausible cases of electron 
distribution function. As the canonical case referred to as 
case (a), we consider two-temperature thermal plasmas, 
where protons and electrons have different temperatures 
and contributions of non-thermal components to the to- 
tal pressure are negligible. As an alternative, we also 
examine case (b) where protons and electrons take the 
same temperature without non-thermal components. 

We further explore two cases (c) and (d) in which non- 
thermal population makes a dominant contribution to 
the total pressure with a negligible pressure of thermal 
population. For the non-thermal population, we assume 
the powe- law distribution functions: 

" NT (7-) oc 7- S<! (7-,min < 7- < 7-, max), 
n p T (lp) <*lp Sp (7p,mm < 7p < 7p,max), (3) 

for case (c) with s p = s e > 2. Observations of the spec- 
tral index in the radio lobe of Cygnus A suggest s e > 2 
(e.g., Carilli et al. 1991; Yaji et al. 2010). 
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Lastly, we set case (d) in which the number spectrum 
of non-thermal electrons is given by a broken power law: 



-S e ,l 



(7-, min < 7- < 7-,crit), 



[7-crit 7- ^ (7-,crit < 7- < 7-,max), 
n p T (7p) « 7~ Sp (7p,min < 7p < 7p,max), (4) 

where s e ,i < 2 and s p > 2 are satisfied. This model is 
based on Stawarz et al. (2007) who suggest that observed 
spectra at the jet termination shock (hot spot) of FR II 
jets (Cygnus A) can be explained by the break at non- 
thermal electron energy (hereafter 7±, cr it)- This type of 
spectra may be due to the absorption of electromagnetic 
waves emitted at the harmonics of cyclotron frequency of 
cold protons, as discussed by Hoshino ct al. (1992) and 
Amato and Arons (2006). Some observations for other 
FR II sources could also be compatible with this picture 
(e.g., Perlman et al. (2010) for 3C445.) 

For cases (c) and (d), the minimum energy of non- 
thermal electrons/positrons (7±,min«ieC 2 ) and protons 
(7p,min?TieC 2 ) are generally assumed as 



7±,: 



(5) 



which is expected when protons and electrons/positrons 
are separately heated and accelerated at termination 
shocks. On the other hand, the values of the max- 
imum energy of non-thermal pairs ('J± t m ax m e c 2 ) and 
protons (7 p ,maxW p c 2 ) are largely uncertain. While 
7±.max™eC 2 may be significantly affected by radiative 
coolings, 7p,maxTO p c 2 may reach the range of highest en- 
ergy cosmic-rays (e.g., Takahara 1990; Rachen and Bier- 
mann 1993). It is reasonable to suppose that 7±,max 3> 

7i,min and 7p,max 7p,min- 

3. TOTAL PRESSURE P 

In this section, we briefly describe the basic idea of esti- 
mating the total pressure P. In Fig. Q]wc show a cartoon 
of the interaction of the jet and ICM. Heating and accel- 
eration processes work at hot spots and those particles 
are injected into cocoons. The cocoon model was pro- 
posed by Begelman and Cioffi (1989) in which the dissi- 
pated energy of jet bulk motion is the origin of the total 
pressure of cocoon and a cocoon of FR lis is expected 
to be overpressured against ICM pressure (Picm) with a 
significant sideways expansion. Therefore, the assump- 
tion of P = Picm is not valid. We have proposed the 
method of dynamical constraint on P by comparison of 
the cocoon model with the actually observed morphology 
of the cocoons (Kino and Kawakatu 2005; 108) and the 
method is applied to various radio lobes (e.g., Machalski 
et al. 2010). We use this model in the present work. 
The reliability of the expanding cocoon model is well ex- 
amined in Kawakatu and Kino (2006). The results of 
relativistic hydrodynamical simulations of Scheck et al. 
(2002) and Perucho and Marti (2007) support the above 
analytic model. The mass and energy injections from the 
jet into the cocoon, which govern the cocoon pressure P 
and mass density p averaged by the source age (i a ge) are 
written as 
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PV 



It, 



2T? 1 A i 



2P 



ft c r>j> ( 6 ) 



P v 



— =2JjA h J: p:l)c. 



(7) 



where j, V, Aj, T® 1 , Jj, pj, and Tj are the adiabatic 
index of the plasma in the cocoon, the volume of the 
cocoon, the cross-sectional area, the total energy flux, 
and rest mass flux, rest mass density, and bulk Lorentz 
factor of the jet, respectively. The term V is evaluated as 
V = 2(ir/3)7Z 2 LS 3 , where LS and 7Z are the linear size of 
the cocoon along the jet axis and the aspect ratio of the 
cocoon, respectively. Here we denote physical quantities 
of the jet with the subscript j. Throughout this work, we 
focus on a relativistic jet. Correspondingly, the shocked 
plasma has relativistic energy, thus we take 7 = 4/3. The 
PdV work done by the cocoon against ICM is taken into 
account in the energy equation Eq. ^ following 108. For 
given picM; we can dynamically estimate total pressures 
P by measuring LS, TZ, and the head cross-sectional area 
of the cocoon. Here the relations of LS — AisCt a gc and 
TZ = l c /LS < 1 hold where l c and /3h s c are the lateral 
size of the cocoon and advance velocity of the hot spot, 
respectively. Since TZ and /?h s have some uncertainties, 
actual P is bounded by maximum and minimum values 

Pmin < P < Pnax- (8) 

Thus we can obtain the total pressure of cocoon P, which 
includes the partial pressures of non-radiating particles. 
The estimate of P has actually been done by 108 for some 
FR II sources and we adopt P values in 108 in this work. 

4. PRESSURE AS A FUNCTION OF r) 

In this section, we express P as a sum of the partial 
pressures and represent it as a function of rj (we call this 
Equation of State, EOS) for respective cases. 

4.1. Case (a) 

First, we examine the canonical case of two- 
temperature thermal plasma. Here we assume that 



pNT _ pNT _ pNT 



and n^ T = n^ T = n^ T = 0. 



The EOS in the cocoon filled with relativistic plasma is 
given by 



F«P ± T 



P T 

p 



= (n T +nl)kT± +nlkT p , 



(9) 



where T±, and T p are the electron/positron tempera- 
ture, and proton temperature, respectively. Hereafter 
we adopt T± = T_ = T+ where T_ and T+ are tempera- 
tures of electrons and positrons, respectively. Following 
Kino et al. (2007), we can obtain T± and T p from Eqs. 
©, 0, and ®: 



kT+ 



Tjin e c 2 



TjiripC 2 



(10) 



which are typically given by kT± =1.3 (jo^j MeV, and 



kT v = 2.3 



GeV. Here we assume the limit of in- 



efficient e/p-coupling i.e, protons and electrons are sep- 
arately thermalized so that kT± = (m e /m p )kT p since 
plasma number densities in large scale jets are conserva- 
tively expected to be too dilute to achieve efficient e/p- 
coupling (e.g., Kino ct al. 2007 and references therein). 
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The emission timescale is so long that radiative cooling 
is negligible. It is worth noting that the geometrical fac- 
tors in Eqs. (|6|) and ([7]) are completely canceled out and 
kT± and kT p are governed only by Tj . 

Inserting Eq. (|10|) into Eq. ([9]), we rewrite the total 
pressure in the cocoon P as 



P(r)) = 2.05 x 10 



- 6 n T 



2 - r? + v— 

m e 



j± ) erg cm 3 



(11) 

where the first term and second term in the square 
bracket correspond to the partial pressure of pairs and 
protons, respectively. 

4.2. Case (b) 

As an opposite extreme to case (a), here we consider 
the case of one-temperature plasma. In this example, 
some of the proton energy is somehow transferred to 
electrons/positrons to achieve an efficient e/p-coupling. 
Then hotter electrons/positrons and colder protons are 
produced. From the condition kT± = kT p , and Eqs. ([6]) 
and 0, we obtain 



kT± = kT„ 



T- 3 m e c 2 



(2 - rj) + 



(12) 



In this case, each population (i.e., p/e~ /e + ) has the same 
kinetic energy. The total pressure is given by Eq. (jlip 
the same as case (a). The essential difference from case 
(a) is that kT± in case (b) is much higher than the one 
in case (a). 

4.3. Case (c) 

For comparison with the canonical case (a), we exam- 
ine case (c) when the cocoon pressure is dominated by 
non-thermal particles. Case (c) concerns when the spec- 
tral indices of non-thermal particle energy distributions 
satisfy s p = s e > 2 as some theoretical work on relativis- 
tic shocks suggests (e.g., Bednarz and Ostrowski 1998; 
Kirk et al. 2000; Achterberg et al. 2001; Spitkovsky 
2008; Sironi and Spitkovsky 2011) and as the radio lobes 
of Cygnus A show s e > 2 (e.g., Carilli et al. 1991; Yaji 
et al. 2010). In this case, electrons and protons with the 
lowest energies are the main carriers of energy. Then, the 
evaluation of partial pressures of non-thermal plasma is 
basically the same as in case (a) when we replace kT± 
with 7±,minTOeC 2 and kT p with 7p j minT- p c 2 . Then P is 
given by 

P(rj) 



r jU^ 1 m e c 2 s e — 1 



(13) 



From this, it is clear that we can appropriately evaluate 
rj for the case (c) by replacing n_ to n^ T in the same 
way as case (a). 

4.4. Case (d) 

Here we examine the pressure of non-thermal electrons 
when they follow a broken power law spectrum Eq. (|4"|). 
Stawarz et al. (2007) indicated 7±. cr it ~ rn p /m e for the 
hot spots in Cygnus A. The energy of the electron com- 
ponent is governed by those with break energy, while 
the number is dominated by those with lowest energies. 



Since s p > 2 is satisfied, lowest-energy protons carry the 
most energy. Therefore, the total pressure P is expressed 
as 



s P i — 1 . , . s r 
I' 1 — A±(2-r,) + ^ 

S e 1 + 2 Sr 



1 n% v 

2 nip 



(14) 



where A± = (7±,crit/7±,min) Se ' 1+2 . Thus r\ can be eval- 
uated when we replace wJL to n^ T and include factor A± . 

5. TESTING PLASMA COMPOSITION 

We explain the method for constraining plasma com- 
position of AGN jets for thermal plasma cases (a) and 
(b) in 5.1, 5.2, and 5.3. The application to non-thermal 
plasma cases (c) and (d) can be readily understood and 
is explained in 5.4. 

5.1. Characteristic pressures 

Firstly we define characteristic pressures which divide 
the number-density /pressure plane into several regions 
as shown in Fig. [2j As a preparation, here we define rj eq 



>k 



m p /m e — 1 



= 1.1 x 10~ 3 (P± = P p ). (15) 



The partial pressure of proton-associated electrons is im- 
plicitly neglected since it is subdominant in the case of 
inefficient e/p-coupling. The line with n_ = 1 x 10 3 n, p 
divides the pair-supported and proton-supported cocoon 
in the limit of inefficient e/p-coupling plasma. By defini- 
tion, the cocoon with ri > n cq is proton-supported (dark 
gray region in Fig. [2]) while the cocoon with r/ < 7y cq is 
pair-supported one (light gray region in Fig. [2]). When 
n_ is bounded by n_. m ; n and n_. max as argued in the 
next subsection, the allowed region of n_ is segmented 
by some characteristic pressures by the characteristic val- 



ues of ri- and rj, i.e., n_ 



c rj = 0, rj = rj CCl , and 



rj = 1. Here, we define six characteristic pressures as 
follows; 



P(r) = 0;n_ = n_ !min ) =P ,, 

P(V = f? cq ;«- = «-,min)=-Pcq 

P(n = 0; n_ = n_ !lnax ) = P ,, 

— ^?cqi ^— ^-,max) — -^eq,i 

P(j] = l;n_ = n_ ;Tnill ) = Pi ;l; 
P(r/ = l;n_ = n_ ?max ) = Pi :lnax . (16) 

Then, by definition, we have the following relations 



,max: 

,max? 
,min j 



-M3,min • -*eq,min 



Po. 



max 



Pe 



;q,max • l,min ■ ^l,max 



Pi 



= 1:2: : 2^^ : ^ : ^V^L j{ll) 

it* — ,min Tl— } niin 7Tl e — 



where we approximate 2 
pressures, we estimate n„, mm 
subsection. 



?7 q ps 2. To evaluate these 
and n_ max in the next 



5.2. Estimation of ri- 

Here we constrain the number density of electrons in 
the cocoon (n_). Wc denote the lower and upper limits 
of ri- as n_.rnin and ti-.n™, respectively. The values of 



,min and Ti- ; max 

show them below. 



are independently constrained and wc 
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5.2.1. Lower limit of n_ 

Here we estimate the lower limit of n_ and examine 
the case when the number density of thermal electrons 
is larger than that of non-thermal electrons nl > ?i^ T , 
since non-thermal electrons are partially injected from 
the background thermal electrons. (Later, the extreme 
cases of < n^ 1 - will also be discussed, being identical 
to cases (c) and (d)). Since the shocked plasma at hot 
spots expands sideways and is injected into the cocoon, 
we can estimate n_ by using where n^ s T is the 
number density of non-thermal electrons in a hot spot. 
We stress that n[J s T is well constrained by observed non- 
thermal emissions of hot spots for FR II sources (e.g., 
Harris and Krawczynski 2006 for review). By connecting 
the number density from the jet to the cocoon based on 
Eq. ([7|) and shock conditions along the jet axis shown in 
Kino and Takahara (2004) ( hereafter KT04), we obtain 



2Vj3 hs 



(18) 



In general, number density of non-thermal electrons with 



power law distribution n^ s T oc J r 
given by 



7hs 
7hi 



7 hs Shsd 7hs can be 



NT -Shs+1 
n hs K 7hs,min • 



(19) 



2 and 7h. 



We assume the standard value of s\ v 

5.2.2. Upper limit of n_ 

The upper limit of n_ can be constrained by the 
absence of thermal brcmsstrahlung from hot elec- 
trons in the cocoon/lobes viewed in X-ray observa- 
tions (Wilson ct al. 2000, 2006). The observed 
X-ray emissions associated with radio lobes arc non- 
thermal emissions and there is no evidence for ther- 
mal X-ray emission from coccons/lobes (Harris and 
Krawczynski 2006 for review). From this, we can 
safely use the condition of £x.obs > L} ilem (n?,T±) 
where L hlcm /V = afr 2 .m e c 3 (n 1 ^) 2 F±(Q±) erg s _1 cm~ 3 , 
F±(0±) = 48e±(hil.ie± + 5/4), and 9± = kT ± /m e c 2 , 
for bremsstrahlung at relativistic temperature (Eq. (22) 
in Svensson 1982) and at and r e are the fine structure 
constant and the classical electron radius, respectively. 
From this, we obtain the maximum n_ as follows: 



Va{r 2 m e c 3 F±(Q±) 



1/2 



(20) 



It is worth commenting on the availability of constrain- 
ing the upper limit of n~ by the analysis of the internal 
depolarization of the radio lobes. Relativistic plasma 
makes a smaller contribution to Faraday rotations since 
electron inertia increases for the relativistic regime and 
it suppresses rotations of polarization angle (e.g., Ichi- 
maru 1973; Melrose et al. 1997; Quataert and Gruzinov 
2000; Huang and Shcherbakov 2011). Therefore, it is not 
effective to use the constraint by RM in the present work. 

5.3. Estimation of n p 



Once n_ is estimated, the proton number density n p 
can be determined as 

n p = r]n-, (21) 

by definition. Here, of course, the conditions of < r] < 1 
and 

^— ,min — W — — — .max 

are imposed. In Fig. [3l the 
allowed region of n p is added to that of n_ shown in Fig. 
[21 In the same way as Fig. [21 the plane is divided into 5 
regions. 

Finally, the allowed regions of n p and n_ can be ob- 
tained by adjoining the range of P. The allowed regions 
drawn in Fig. [3] are bounded by Eq. ([8]). Thus, we can 
obtain the definitive allowed regions of n p and n_. 

5.4. Application to cases (c) and (d) 

In 5.1, 5.2 and 5.3, we consider physical quantities as- 
sociated with thermal plasma in cases (a) and (b). But 
those can also be applied to non-thermal plasma by the 
proper replacements of number densities and average en- 
ergies of particles. With regard to average energies, we 
have already explained the replacements in the previous 
section. As for n_, m i n , the estimate shown in 5.2.1. can 
be applied both for thermal and non-thermal plasma. As 
for ft-, m ax, the estimate shown in 5.2.2. can be applied 
only for thermal plasma. So, we do not use n_ jmax for 
the cases (c) and (d). Thus we can properly estimate rj 
for cases (c) and (d). 

6. APPLICATION TO CYGNUS A 

Here we apply the above method to Cygnus A (z = 
0.0562) which is one of the best studied FR II radio 
galaxies (e.g., Carilli and Barthel 1996; Steenbrrugge et 
al. 2008, 2010; Yaji et al. 2010). The physical quanti- 
ties of Cygnus A have been well constrained by previous 
work. To constrain the real values of P and n_, we care- 
fully evaluate TZ, /3hs, and Tj. The term TZ has an effect 
on n_ via a cocoon volume V. The term Tj is directly 
proportional to P. The term /3h s controls the source age 
i agc which governs the injection rates of mass and en- 
ergy into the cocoon. These arc summarized in 6.1. The 
resultant allowed region of n_ and n p is summarized in 
6.2. 

6.1. Viable ranges of physical quantities 

We show adopted conditions of the model parameters 
for deriving the above results. We fix the cross section 
area of the jet as Aj = nR 2 s = tt(2 kpc) 2 (Wilson et al. 
2000) and the number density of ICM just ahead of the 
hot spot as tiicm = 0.5 x 10~ 2 cm" 3 (the shell No. 6 in 
Table 5 in Smith ct al. 2002). 

• Cocoon morphology TZ. 

From images of the Cygnus A cocoon, we can di- 
rectly constrain TZ. The upper limit 1Z ~ 0.5 is de- 
termined by the Chandra X-ray image (Wilson et 
al. 2000, 2006; Yaji et al. 2010). The lower limit 
TZ « 0.25 is directly measured by the 330 MHz 
VLA image (see also Carilli et al. 1991; Lazio et 
al. 2006). Therefore we set 

0.25 < TZ < 0.5, 

in the present work. 
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Cocoon head velocity /3hs- 

Cocoon head velocity which equals the hot spot 
advance velocity (/?hs) is well constrained by the 
synchrotron aging method. The estimated /3hs 
has some uncertainty due to the uncertainty of 
magnetic-field strength in the cocoon. From the 
result of synchrotron aging diagnosis in Carilli et 
al. (1991), we adopt the allowed range of /3hs as 

0.01 < /3 hs < 0.06. 

We emphasize that sufficiently large uncertainty is 
taken into account here. The adopted value of /3hs 
is quite typical for hot spots in FR II radio galaxies 
(e. g., Scheuer 1995). 

Lorentz factor of the jet Tj . 

It is difficult to determine the true velocity of the 
jet. At least we may say that apparent velocity of 
blobs obtained by VLBI observations show a min- 
imum velocity of underlying flow. A fast apparent 
motion of a blob at the jet base (0.56 ± 0.28) c has 
been reported by VLBI observations (Bach et al. 

2003) . Furthermore, suggestions of superluminal 
motion were made (Krichbaum ct al. 1998; Bach 
et al. 2002) although they had not been clearly 
confirmed. On VLA scale, a clear asymmetry in 
brightness distribution of a kpc-scale jet due to 
a relativistic motion is seen (Perley et al. 1984). 
Therefore, overall radio observations seem to in- 
dicate relativistic motion. Bearing this in mind, 
we assume that the jet is relativistic and the four- 
velocity of the jet Fj/Jj is set as 

i < r^j < 30. 

Here the upper limit is assumed as Tj « 30 based 
on the statistical study of radio jets of MOJAVE 
sources (Lister et al. 2001, 2009; Kellermann et al. 

2004) . 

Cocoon pressure P. 

Using the value of V = 1 X 10 7Q 1Z 2 cm 3 , we can 
estimate the total pressure P as 

x 10" 11 erg enr 3 < P < 4 x 10~ 9 erg cm" 3 . (22) 

The lower limit equals the ICM pressure 8 x 
10" 11 erg cm -3 measured by Arnaud et al. (1984) 
to satisfy the over-pressured cocoon condition. Al- 
though the upper limit of P is basically adopted 
from 108, the value 4 x 10~ 9 erg cm" 3 is 4 times 
larger than the original estimate in 108. This is 
due to the change in minimum value of 1Z from 
0.5 to 0.25 based on VLA's 0.3 GHz image. It 
should be stressed that our adoption of the allowed 
range of P is sufficiently wide compared with all 
of the previous work (e.g, Carilli 1998 for review). 
Note that Yaji et al. (2010) estimates that P^ T in 
the radio lobes as P^ T « (1 — 2) X 10~ 9 erg cm" 3 
for 7± ~ 1 which causes P^ T > P m i n . So, if P 
completely equals the radio lobe pressure, then the 
range P m i n < P < P^ T is excluded and the allowed 
P range becomes narrower. The allowed example 
with P mm < pNT < p < P max is involved in cases 
(c) and (d). 



• Non-thermal electron number density n_ hs . 

The lower limit ?i_. m i n largely depends on n^ s T . For 
Shs = 2, the number density of non-thermal elec- 
trons in the hot spot can be obtained from 

n^ixio-^^)" 1 cm -3 ; (23) 

via detailed comparisons of the SSC model with the 
observed broadband spectrum (Wilson et al. 2000; 
KT04; Stawarz et al. 2007) where 7 hs ,min « Ij. 
We stress that these three independent papers de- 
rive similar values of n^ s T although Stawarz et 
al. 2007) adopts the different electron-distribution 
function shown in Eq. Furthermore, we 

note the importance of low-frequency radio spec- 
tra since it affects the estimate of n^ s . Regarding 
low-frequency radio observation, we briefly com- 
ment on the work of Lazio ct al. (2006). They 
indicated spectral flattening and turnover at ~ 
100 MHz. However it seems difficult to deter- 
mine these accurately because the spot sizes are 
smaller than the VLA beam sizes at the above fre- 
quencies. The LOw Fr equency ARray (LOFAR) 
(http://ww w.lofar.Org/l) an d Square Kilom eter Ar- 
ray (SKA) (http: / /www.sk atelescope.org/[ ) will, in 
future, tell us the real turnover frequency with suf- 
ficiently high resolution. 

• Thermal electron number density . 

Here we comment on the difficulty of constraining 
n^_. We use the absence of brcmsstrahlung emis- 
sion. The A-ray observations for Cygnus A show 
the flux upper limit as ~ 1 x 10 -13 erg s" 1 cm" 2 
(e.g., Smith et al. 2002). 

As already mentioned, the constraint from the in- 
trinsic RM is not available, because plasma tem- 
perature is relativistic in the present work. Even 
worse, Cygnus A is known for its unusually large 
RM values and thus it is not a good example from 
which to argue the intrinsic depolarization (Dreher 
ct al. 1987; Garrington and Conway 1991). No 
evidence for intrinsic depolarization between 5 and 
15 GHz is found and the origin of the large RM is 
thought to be the external bow shock which sur- 
rounds the radio lobes (Dreher et al. 1987; Carilli 
ct al. 1988). Hence it is not appropriate to use the 
constraint from RM for Cygnus A. 

6.2. Results 

Below we show resultant allowed region of tt,_ and n p 
for cases (a), (b), (c) and (d). 

6.2.1. Case (a) 

Considering the uncertainties of Tj/Jj and /3h s , we exam- 
ine two limiting cases with Tj/^j = 1 and /?hs = 0.01 being 
a High-n case, and that with Tj/3j = 30 and /3hs = 0.06 
being a Low-n case. For the High-n case, n_ is about 
two orders of magnitude larger than that of the Low-n 
case. 

In Fig. 21 we show the allowed region of n_ and n p 
for the High-n case. First of all, we find that n_ > n p 
always holds and this satisfies ry~10" 2 atP = Pmax- 
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This implies that positron mixture is inevitable. In other 
words, Pi^iin is much larger than P max obtained by the 
Cygnus A cocoon calorimetry. (If we are force to make 
-Pi,min smaller, then 7 m ; n becomes larger and such a case 
coincides with (b).) The allowed regions of n_ and n p 
arc further divided by two regions. The pair of light- 
gray regions show the one in which P± > P p is satisfied. 
On the contrary, the pair of dark-gray regions display 
the one in which P± < P p holds. Interestingly, we find 
that the regions of P p < P± and P p > P± are both wide 
in the range of allowed P. Only in the range of P ~ 
(3 — 6) x 10 -10 ergem -3 , the pair dominance P p < P± 
alone is permitted in the High-n case. 

Fig. [5] displays the result for the Low-n case. Simi- 
lar to the High-n case, n_ > n p always holds and they 
satisfy r\ ~ 10 _1 at P = P max . Due to the decrease 
in n_. m ; n , the number densities in allowed regions are 
about two orders of magnitude smaller than that for the 
High-n case shown in Fig. |U Correspondingly, Po, m in, 
Peq,min: and -Pi^min decrease. Since P e q,max < Pnax is 
still satisfied, both of the regions with P p < P± and 
that with P p > P± are allowed in this case. In other 
words, the Low-n case also draws the same conclusion 
with the High-n case qualitatively. Quantitatively, the 
upper limit of n p becomes larger when n_ jm ; n becomes 
smaller and correspondingly the maximum n achieved at 
Pnax becomes larger by a factor of ~ 10 than that for 
the High-n case. 

In summarizing case (a), we find that n < 1 always 
holds in the allowed range of P. In other words, this 
indicates the existence of pairs in the cocoon. We 
find that (i) e^- pair is dominant in terms of number 
density, and (ii) both the "pair-supported cocoon (i.e., 
P± > Pp)" and the "proton-supported one (i.e, P± < 
P p )" are allowed. The pair-supported cocoon is different 
from the previously suggested one in which protons are 
dynamically dominated (e.g., De Young 2006). 

6.2.2. Case (b) 

For Cygnus A, we face a difficulty of realizing one- 
temperature plasma. First, let us consider the case of 
same n_. m i n as in Figs. OandUl Then all of these thermal 
electrons should be heated up to kT± ~ 10 4 m e c 2 and 
injected into the lobes in the case (b). In the radio lobes, 
Yaji et al. (2010) evaluates the number density of non- 
thermal electrons as ~ 10~ 7 cm" 3 at 7_ ~ 10 4 . So, if 
we allow the existence of thermal plasma with the same 
"■-.min in Figs. [4] and [5] but with kT± = kT p ~ 10 m e c 2 , 
a big thermal bump at ~ 10 9 Hz should appear. However 
there is no such bump in the observed spectra of the 
radio lobes. Therefore, we can exclude the case of the 
same n_ jm j n with kT± = kT p ~ 10 4 m e c 2 . 

Next, we consider smaller n_ m i n . Using the relation 

rc.-, m in °c 7i7slmin in E< 1- 033), tnc increase in 7 hs , m in leads 
to the decrease in n_ jln i n in Figs. |4] and basically, 
7hs,min ~ 10 4 is required at the hot spot (e.g., Harris 
et al. 2000; Hardcastle, Birkinshaw, and Worral 2001; 
Blundell et al. 2006; Godfrey et al. 2009). However, 
in the case of Cygnus A, the model spectra of the hot 
spots with 7hs,min > 2000 conflict with the observed ones 
(KT04). Therefore, case (b) is not likely for Cygnus A. 

6.2.3. Case (c) 



Let us consider the case of dominant non-thermal pres- 
sures and a separate acceleration of electrons and protons 
with a steep power law spectrum. This is almost iden- 
tical to (a). A slight difference between this case and 
(a) is the evaluation of n_. m ; n . Since non-thermal pairs 
are dominated in this case, the allowed region would be 
limited around n_ m n_. m j n in Figs. 2]and[5j 

6.2.4. Case (d) 

Let us consider case (d). The factor A± = 
(7 ± crit /7 ± min )~ s + 2 in Eq. (fLfj) is the only element to 
change the result from (a). Since 7 cr it,± ~ m p /m e is sug- 
gested by Stawarz et al. (2007), we can estimate A± as 
A± « 14(rj/10) - 5 for s e i = 1.5. Therefore, a difference 
between this case and (a) is the larger P± by a factor of 
A± . Although the spectral break may be suggested from 
radio observations for case (d) , n^ T is dominated by elec- 
trons at a break energy 7 cr it,± TO eC 2 and proton energies 
are not entirely transported to electrons. Therefore, re- 
sults of (d) are expected to be intermediate between cases 
(a) and (b). 

7. SUMMARY AND DISCUSSIONS 

In this work, we propose a new method for testing 
plasma composition of AGN jets by using the cocoon 
dynamics. In particular, we properly evaluate partial 
pressures of protons and pairs. The point of the 
method is that n p and P p can be constrained by con- 
sidering the global conservations of kinetic energy, mass, 
and momentum of shocked plasma in the cocoon. Re- 
garding particle distribution functions in the cocoon, it 
is hard to determine them uniquely because of sparse- 
ness of observational data. Therefore, we examine four 
typical cases in this work. Cases (a), (b), (c) and (d) re- 
spectively present two-temperature thermal plasma, one- 
temperature thermal plasma, non-thermal plasma with 
their spectral indices harder than two, and non-thermal 
plasma with a broken power-law electron spectrum. 

The three significant advantages of the present work 
compared with previous work are summarized as follows; 

1. P estimate is based on global cocoon dynamics. 
Since it is beaming-independent calorimetry of the 
true amount of energy released by the jet, the esti- 
mate of P from cocoon dynamics has fewer uncer- 
tainties compared with blazar studies. 

2. We focus on powerful FR II sources. Rclativis- 
tic hydrodynamic simulations tell us that FR II 
sources have less entrainment phenomena than FR 
I sources. Therefore, FR lis are better for testing 
genuine plasma composition of AGN jets. 

3. We properly deal with the partial pressure of ther- 
mal electrons / positrons P± . Although P± is a crit- 
ically important finite quantity, most prior efforts 
assume PjF = merely for simplicity. 

Applying the method to the best studied FR II source 
Cygnus A, we draw the following conclusions which pri- 
marily indicate the existence of numerous e ± pairs in the 
cocoon of Cygnus A. 

• Cases (a), (c) and (d), in which the average energy 
of electrons and positrons is significantly lower than 
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that of protons (77 < 10 _1 for Low-n case; r] < 10~ 2 
for High-n case), are allowed without violating the 
observational constraints. The results in (a) and 
(c) are almost the same, except that the lowest 
energy electrons are thermal ones and non-thermal 
ones for (a) and (c), respectively. Cases (a) and (d) 
also show similar results but for a larger P± in (d) 
by a factor of ~ 14 than the one in (a) . 

• We can rule out case (b) in which electrons and 
positrons are heated up to the proton temperature 
of ~ 10 4 m p c 2 . Because there is no thermal bump 
due to the hot thermal plasma. 

• For (a), (c) and (d), we find that the number 
density of is larger than n p in any allowed P 
and the obtained n + is always more than 10 times 
larger than n p . We conclude that pure e/p plasma 
is excluded and e -proton mixture composition is 
achieved in the Cygnus A jet. Therefore, further 
studies on the e^ pair loading problem extending 
previous ones (e. g., Blandford & Levinson 1995; 
Li & Liang 1996; Thompson 1997; Beloborodov 
1999; Yamasaki, Takahara & Kusunose 1999) will 
be more important and the study of its bulk accel- 
eration of e outflow (Iwamoto and Takahara 2002, 
2004; Asano and Takahara 2007, 2009) will also be 



highly motivated. 

• We find that both e/p plasma and e^ pair pressure 
supported scenarios are permitted within the limit 
of current observational constraints. We quantita- 
tively show the allowed regions of P p > P± and 
Pp < P± by our new method (see Figs. H]and[5]). 

Lastly we add a brief comment on Pp T . Recently 
Atoyan and Dermer (2008) has suggested the possibility 
of a secondary emission induced by high-energy protons 
at Cygnus A. The luminosity of the secondary emission 
depends on Pp T . If the emission is detected in the fu- 
ture, it will provide us a new direct constraint on Pp T . It 
could also give us a new constraint on cosmic-ray propa- 
gations influenced by the galactic magnetic field (Dermer 
et al. 2009). 
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Fig. 1. — A cartoon of a powerful FR II radio galaxy. A pair of jets is ejected from the core and they are decelerated via strong shocks. 
The shocks are identified as the hot spots and the remnant of decelerated jets envelopes the overall jet system and this is identified as a 
cocoon. Part of cocoon is normally observed as radio lobes. The cocoon head and the hot spots advance at a speed v-^g. Swept-up ambient 
matter becomes a shell and surrounds the cocoon. The projected linear size is denoted as LS in this work. 
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Fig. 2. — A schematic diagram of the allowed region of n_(P) plotted versus the cocoon pressure P for given Tj. They are limited 
by n — ,min !^ n - < "-.max and < < 1. The region of e^-supported cocoon (P± > P p ) is colored in light gray while the region of 
proton-supported cocoon (P p > P±) is colored in dark gray. 
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Fig. 3. — Same as Fig. [2] but including the allowed region of n p . When rj = 1, n_ = n p holds by definition. When r) < 1, positron 
mixture is required by the charge neutrality condition of n_ = n p + n+. The plane is divided into 5 regions by characteristic pressures. 
The actual allowed region is further limited within P mln < P < Pmax by the consideration of cocoon calorimetry. 



B 




P (erg cm ) 

Fig. 4. — Allowed regions of n_ and n p for Cygnus A with Tj/3j = 1 and /?hs = 0.01 (we call this the High-ra case). The region within 
8 X 10 — 11 erg cm~ 3 < P < 4 X 10 — 9 erg cm~ 3 shown here is the one allowed for Cygnus A. As explained in Fig. [3] the region in which 
P± > P p holds is colored in light gray while the region where P± > P p is satisfied is colored in dark gray. It is found that er 1 pairs always 
dominate in terms of number density but either "pair-supported cocoon (i.e., P± > P p )" or "proton-supported one (i.e, P± < P p )" is 
possible. 
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Fig. 5. — Same as Fig. [5]but with Tj = 30 and f}^ = 0.06 (we call this the Low-n case). Although the allowed regions of n_ and 
are about two orders of magnitude smaller than the ones in Fig. [4] (High-n case), the Low-n case also draws the same conclusion as 
High-n case. 



